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Abstract 

We study the Poschl- Teller equation in complex domain and deduce infinite families 
of TQ and Bethe ansatz equations, classified by four integers. In all these models 
the form of T is very simple, while Q can be explicitly written in terms of the Heun 
function. At particular values there is a interesting interpretation in terms of finite 
lattice spin--^^ XXZ quantum chain with A = cos J- (for free- free boundary con- 
ditions), or A = — cos Y (for periodic boundary conditions). This result generalises 
the findings of Fridkin, Stroganov and Zagier. We also discuss the continuous (field 
theory) limit of these systems in view of the so-called ODE/IM correspondence. 
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1 Introduction 



Some years ago, an unexpected connection was found between certain + 1 dimensional 
quantum-mechanical problems and 1 + 1 dimensional conformal field theories j1j|[2t|31|^ 
ISlinilZllH] • The simplest example involves the Schrodinger equation 

--^H^, E) + (x" - E)i;{x, E)=0, (1.1) 

and the fact that ()1.1|) has unique solution y{x,E), entire in x and E, which decays 
along the positive real axis as x ^ oo. The function y(x, E) can be shown ^ to satisfy 
a Stokes relation 

T{E)y{x, E) = uj-^l'^ y{ujx, uj^'^E) + w^/^ y{uj'^x, co'^E) , (1.2) 

with UJ = exp(z27r/(a + 2)). This is strikingly similar to the TQ relation, a functional 
equation which was introduced in the context of the six- vertex model by Baxter (cf . ^ ) . 

The correspondence is actually much more precise: the Stokes multiplier ^U] T{E) 
and the spectral determinants Q^{E) = y{x,E)\x=o and Q~{E) = y'{x,E)\x=Q are 
equal to the ground-state eigenvalues of the transfer matrix and of the Q operators, 
respectively, of the six-vertex or, equivalently, of the spin-^ XXZ quantum chain, in 
suitable continuum (field theory) limits. The continuum limit of the six-vertex model 
is related to conformal theory at c = 1, and the description of such theories in a 
framework similar to Baxter's lattice setup was addressed in an important series of 
works by Bazhanov, Lukyanov and Zamolodchikov ,12, 13. 14^. The correspondence also 
applies, directly, in this context. Since the initial observation of PP, many mathematical 
aspects of the correspondence have been clarified [HISIIS]) but the physical reasons for 
its existence remain mysterious. In addition, up to now, the correspondence has been 
limited to the ground-state energy of the conformal field theory and it is unclear whether 
it can be extended to massive theories, to excited states*, or to finite lattice systems^. 

The purpose of this paper is to report some progress on the last of these questions. 
The key difference in finite lattice problems is the appearance of an extra function ^{E) 
in the TQ relation: 

T{E)Q{E) = b-^ ^{ujE)Q{uj~^E) + h ^{uj'^E)Q{uj'^E) . (1.3) 

(6 is a pure phase.) This function is fixed by the problem under consideration. In par- 
ticular, it encodes the number of lattice sites, and tends to 1 when a suitable continuum 
limit is taken. 

'While we were finishing writing this article the preprint 1151 appeared. In this very interesting 
paper, Schrodinger equations corresponding to excited states are proposed. 

^An extension of the finite lattice Baxter TQ relation with extra coordinates like parameters has 
been recently introduced by Weston and Rossi in for q generic (see also |17I18| for q a root of unity). 
The relationship between these results and the ODE/IM correspondence still needs to be clarified. 
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We start with the fact that at the points a = 2 and a = 1 the solutions of can 
be given exphcitly, in terms of hyper geometric functions. In the search for finite-lattice 
generalisations of the correspondence it seems natural to begin with these two models. 
At the linear point, a = 1, there is the additional advantage that the roles of E and x 
are interchangeable: setting z = x — E, (|1.H1 becomes the standard Airy equation 

-^Hz)+zil:{z)=0 , (1.4) 

and therefore jHQ] 

Q+(^)U=i = Ai(-^) and Q-{E)\^=i = -Ai'{-E) , (1.5) 

which means that the functions Q^{E) themselves satisfy differential equations. 

The ideal situation would be that the finite-lattice version of the problem would 
share this {E, x)-democracy. Surprisingly, this was precisely the discovery made by 
Fridkin, Stroganov and Zagier in |19j : though they did not make a connection with the 
earlier results of ^HHSHHISl , they empirically discovered that the Q function for the 
spin-^ XXZ quantum chain with free- free boundary conditions and A = 1/2 was related 
to a specialisation of the Poschl- Teller equation 

(P 9n(2n + l) , , , , , , , 

^ ' ' X{s) = -X{s) ■ (1.6) 



2cosh2 3s/2 

In a subsequent paper j^H], the same spin chain, but with A = —1/2 and with periodic 
boundary conditions, was also related to the Poschl- Teller equation with a different 
eigenvalue. 

In this paper we show that the most general Poschl- Teller equation (given in equa- 
tion ()2.1() below) contains infinite families of finite TQ and Bethe ansatz equations, 
selected by fixing the four parameters (M, N, L, m) to positive integer values. At the 
particular values (0, L, 0) and (0, A'^, L, 1) there is a straightforward and interesting 
interpretation in terms of finite lattice spin--^^^ XXZ quantum chains. 

The plan of the paper is as follows. The relevant equation and its analytic properties 
are discussed in section |2j The TQ relation is derived in section |21 and the relationship 
with the quantum spin chains discussed in section^ In section [5] some numerical results 
are reported and in section El the continuous limit of the equation is studied in view of 
the ODE/IM correspondence. Finally section [7| contains our conclusions. There are two 
appendices: in appendix the solution in terms of the Heun function is derived, while 
the locations of the trivial zeroes of the solution are discussed in appendix IbI 



2 The differential equation 

We consider the generalised Poschl- Teller equation 

/ N(N + 1) M(M + 1)\ , , y , , , , 
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As explained in appendix^ the differential equation can be mapped into the Heun 
equation, allowing any solution of H2.1() to be written in terms of the Heun function H. 
This is one of the reasons why the Poschl- Teller equation has historically played an 
important role in the quantum-mechanical modelling of two-body problems with short- 
range interactions. In these applications the wavefunction is usually defined on the real 
axis, and the physical requirement of square integrability constrains a to integer values, 
allowing the solution to be written in terms of the more standard 2^1 hypergeometric 
function. 

In this paper we shall instead consider equation (|2.1j) on the whole complex plane, 
and one of the requirements placed on its solutions will be meromorphicity. The demand 
that x{s) be single- valued around the singularities of cosh~^ s and sinh~^ s immediately 
restricts the parameters N and M to integer values. However, in the following we shall 
impose further conditions, which emerge as follows. 

Introduce a new variable = — exp(2s) , and set 

m = aL-l, i;M,N,m{^) = x^"'+^^/\{ln{^/^)) . (2.2) 

Then 'ipM,N,m,{x) is solution of 

/ 2 d L^N{N + l)x^ L^M{M + l)x^ 

"^""di {x^ - 1)2 + (x^ + 1)2 

The requirement that ipM,N,mix) be single- valued on the whole complex plane leads 
to the quantisation of the four parameters A^, M, L, m to integer values. To see the 
quantisation of m, notice that the points x = and x = 00 are, in general, singular 
points of 1)2. 3|) . and in their vicinity solutions behave as 

V'M.TV.mlrc) ~a + /?x™+i + ... , (2.4) 

and therefore single- valuedness constrains m € Z. The case m > is already very 
rich in structure, and so we shall restrict ourselves to this case. Without any further 
loss of generality, we conventionally set L,N,M > 0. All of this was to allow the 
single- valuedness and hence meromorphicity of the solutions to 1)2. 3j) . To single out 
one particular solution, to play the role of the function y in (jl.ip . we shall impose the 
boundary condition 

V'M,iV,m(x)U^i ~ (l-x)^+^ (2.5) 
This condition is natural, in that 1)2. 3|) has regular singularities at 

(x^Tl)U=x,,±=0. (2.6) 

Near these points a generic solution behaves as 

'4>M,N.m{x) ~ c(x - Xj_)"^ , {X^ ~ 1) , (2.7) 
V'M,iV,m(x) ~ C(X - Xi+)-^ , (X^ ~ -1) . (2.8) 



j 1pM,N,m{x) = . (2.3) 
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If we impose c = then, exceptionally, 



i^M,NAx) ~ d{x - , (2.9) 

which up to the (arbitrary) normalisation is exactly the condition (|2.5j) . The choice 
to impose boundary conditions near regular singularities in this way might seem to be 
unmotivated at this stage, but it will be crucial in making a connection with in 
the scaling limit. We shall return to this point in section |H1 

Notice that equation (|2.3() is invariant under the transformation 

(x, M, N) {xuj^/'^.N, M) , (2.10) 

where ui = exp(2i7r/L), and consequently also under 

(x, M, N) {xio, M, N) . (2.11) 

These symmetries force further zeroes in tpM,N,m{x). Being images of x = 1 under 
certain rotations, they are located on the unit circle, and are, in some sense, trivial. 
They will, however, contribute non-trivially to the Bethe ansatz equations which fix the 
nontrivial zeroes - see, for example, 1)3. 11() below. The determination of the locations 
of the trivial zeroes is simple but technical, and we relegate it to appendix IbI 



3 The connection formula and the TQ relation 

We shall now formulate the problem in a setup similar to that used in presence of Stokes 
sectors ^Q]- Set g{x) = i^M,N,m{x) and define 

gk{x) = g{iu''x) . (3.1) 

Then the symmetry 1)2. 11() ensures that gi{x) and g^i{x) are also solutions of equation 
(|2.3|) . Near x = 1 they behave as 

5i(x) ~c+(x-l)-^ , g_i(x) ~c_(x-l)~^ , (3.2) 

and so the pair of functions {go{x), gi{x)} is (apart from the particular values of m = 
L — 1 modL) a basis of solutions. Expanding g^i in this basis and rearranging, 

W[-l,l]go{x) = W[-l,0]gi{x) +W[0,l]g-i{x) , (3.3) 

where the Wronskian is 

W[i,j] = gi{x)g'j{x) - g\{x)gj{x) . (3.4) 

Because of the term — rnx '^^^j^'* in ()2.3() one can deduce that the Wronskian between any 
pair of solutions has the form 

W[g,f]=cstx^ , (3.5) 
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and one can factorise x™ out of (|3.3|1 . We can now use the large x asymptotic (j2.4l) 

g{x) ~ a + 6x'^+^ + . . . (3.6) 

to find the exact expression for 1, 0] and VF[0, 1]. For a, 6 7^ the result is 

W[-l,Q] = 2m6(m + l)a;(™+^)/2 gij^(m+i^)^m ^ ^3^7^ 

T^[0,1] = 2m6(m + l)u;-(™+i)/2 vr)x'^ , (3.8) 

= 4M6(m + l)sin(^7r)cos(^7r)x'" , (3.9) 

and H3.3|) becomes 

m+l m+1 

2cos(^7r)5(x) = w 2 c,(u;"^x)+w~ 2 (/(wx) . (3.10) 

This equation is almost identical to Baxter's TQ relation, save for the fact that Q{x) 
is, by definition, entire while (j{x) = ipM,N,m 

(x) is not. This can be simply overcome by 

introducing a new function q{x) defined as 

e-i Nf, 
(x^ - l)^(x^ + l)''g{x) = \{{x - {Jyf''+^ J{{x - {u:)"-^f'^'+\{x) . (3.11) 

j=0 i=l 

In the above we have used the knowledge of the trivial zeroes and poles discussed in 
appendix El Note that ^, w', ki and are respectively defined in ()B.5|) . (|B.6|) . ()B.7|) 
and at the very end of appendix IbI 

By a consideration of the possible singularities and trivial zeroes in the previous 
section and in appendix El we immediately deduce the following factorised form for 
q{x): 

q{x) = n (1 - ^) ' (3-12) 

i=i ^ 

with 

Xj = l/xK+i-j ■ (3.13) 

The function q(x) also satisfies a TQ-type relation, and it is entire. 

The number of nontrivial zeroes, K, of q{x) is easily evaluated by noticing that 
''pM,N,m{x) is a meromorphic function of x. Then, the asymptotic behaviour (|3.6|) indi- 
cates that Max(m + 1,0) should be equal to the total number of zeroes minus the total 
number of poles existing at finite x. For m + 1 > 0, this leads to 

K = {m + l)+N{L-i) + M{L - Nk) - 1{N + 1) - Nk{M + 1) . (3.14) 

The set of numbers {xj} constitutes the nontrivial zeroes of the wavefunction. They 
are fixed by the Bethe ansatz equations. To match the standard notation, we change 
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variables x — > —x and Ej = —xj, and also set x = exp(u), Ej = exp(Mj). In these new 
variables the connection (TQ) formula p.lUf) becomes 

r(n)Q(n) = (t){u - 2irj)Q{u - 2irj) + (j){u + 2irj)Q{u + 2ir]) , (3.15) 

where we defined 

K 

Q{u) = J]sinh(^^), (3.16) 

j=0 i=l 

t{u) = (-l)^+*^2cos(^7r)0(u) , (3.18) 

with ri = tt/L and r/' = ir/i. 

Prom its explicit form, t{u) should be pole- free while the formal solution of the 
above algebraic equation seems to possess poles at the zeroes of Q{u). Thus the residue 
at these points must be vanishing. This is exactly the same reasoning which leads to the 
Bethe ansatz equation in integrable systems. A suitably-chosen solution to the resulting 
Bethe ansatz equation characterises Q{u), and consequently q{x), and exhibits several 
interesting patterns depending on the choice of parameters. Before presenting examples, 
however, we shall discuss the connection of our findings to quantum magnets. 



4 Model identification 

Consider a one-dimensional quantum system in which quantum spins of magnitude S 
are assigned to each site of a length J\fs chain. They interact via spins of magnitude 1 /2 
living on bonds. The strength of the interaction is characterised by A = cos A. Assume 
further either periodic boundary conditions (p.b.c), or free- free (f-f) boundary conditions 
with Uq{sl2) invariant interaction and q = expzA. For example, the Hamiltonian for 

5 = ^ , and periodic boundary conditions, is: 

n=l ^ 

The transfer matrix Ts{u) is given by either a single (p.b.c) or doubled (f-f) form via 
Skylanin's construction The auxiliary space has spin ^ , and the quantum space is 
given by the A/5-fold tensor product of a spin S space. Then the following TQ relation 



6 



holds: 



T^''\u)Qsiu) = (l)siu-2iSX)Qsiu + 2iX)+(l)s{u + 2iSX)Qs{u-2iX) , (4.1) 

frijsmh^ , (p.b.c) 
I n, sinh ^ sinh (f-f) 



Qs{u) 

4>s{u) 



'nfjisinhi^, (p.b.c) 
^ smh{u) Uaii sinh sinh , (f-f) 



where oja stands for the inhomogeneity, and (u) stands for the renormalised transfer 
matrix: 

Tir)(^^ ^ l^siu) (p.b.c) 
^ I sinh nT5('u) (f-f) 

For periodic boundary conditions, the above relation can be shown directly, while 
for free- free boundaries, it generalises established results for S = ^ and 1 |22ll23j . 

The similarity between ()4.1() and our connection formula p.lSMTTH)) is clear. To 
check the precise correspondence, we now examine some simple examples, taking Af = 
and m = or 1. For m = 1, we additionally impose that L be odd, so that in all cases 
i = 1. Then (p in (||-i.l7j) simplifies considerably: 



4>{u) 



Noting also the property 



cosh2^+i(f) , iNk = 0) 

isinh(u)cosh2^(f) , (Nk = 1) . 



Uj u u 

cosh(- ± if) = ±i sinh(- T i^TT) = ±i sinh(- ^ i{L - 2)2) , (4.2) 

for Nk = 1 it is immediately seen that the connection rule ()3.15M^?TH)) coincides with 
(|4.1jl for the spin chain with f-f boundaries, an even number of sites A/5 = 2N, 
and with parameters r] = X = j^, Uj = Vj, and u)a = 0. For the match to be complete 
the function t{u) should be related to an eigenvalue T^\u) of the transfer matrix as 
2(-l)^+V(u)=r('-)(n). 



Similarly, for A^^ = and L odd, the connection rule coincides with (|4.1j) for the spin 
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^ ^ chain, but with p.b.c and an odd number of sites Ms = 2N + 1. The parameters 
must be identified as A = vr — 77 = -^^^77^, Uj = vj + vri, (— 1) 2"+™~^r(u) = Tg^\u + ni 
and LOa = 0. 

Specialising to L = 3, the above results recover the findings of jl9ll2Uj . The co- 
incidence between the ODE and the spin chain was checked numerically for A^^. = 0, 
L = 3,5, with = 1,2,3. We adopted a 'brute force' diagonalisation of the transfer 
matrices associated to the spin systems, and then verified that the resultant spectra 
contain eigenvalues of the form T^\u) = 2(— l)^+^r(n). These correspond to the par- 
ticular solutions of the Bethe ansatz equations 1)5. If) which will be associated with the 
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ODE in the next section. The eigenvalues are not particular members of the spectra: 
they are neither the largest in magnitude nor the smallest. However, the same BAE 
patterns do play a distinguished role in a particular (isotropic) fused model. See the 
discussion in section [Tj The remarkable simplicity of the expression for Ts{u) comes 
from the elementary expression for t{u) in (|3.18|) : it reflects the special nature of the 
points we are examining even before the scaling limit is taken. 

When Nk = and L is even, our connection rule differs from the periodic boundary 
condition case of 1)4. 1() by a sign. This suggests the need for a different choice of boundary 
conditions for the spin model. We leave this for future work, as well as the identification 
of the connection rule and the TQ relation in higher spin chains, with general choice of 
N, M, L, m, where the inhomogeneities uja should be chosen properly. 

We make one further, more general, remark in closing this section. There are some 
ambiguities in the choice of functions in the lattice model and the ODE : in particular, 
equation is invariant if r(n) and </>(u) are multiplied by a common arbitrary entire 

function of ^. Although the choice we have made above appears to be the most 
natural, and is supported by our numerical results, we cannot exclude the possibility 
of the extra factor being relevant in a more general situation. We hope to be able to 
resolve this issue in a future publication. 



5 The Bethe ansatz equations and string like solutions 

From equation 1)3. 15() . and the reasoning given after that equation, the zeroes {uj} of 
Q{u) satisfy the following Bethe Ansatz equations (BAE)^ 

<P{n,-2^r,)Q{n,-2^r,) _ ^ ^^^.^...^j^y (5.1) 



4>{uj + 2iri)Q{uj + 2iri) 

The solutions to these equations which are related to our ODE exhibit various interesting 
patterns of zeroes depending on the choice of {M, N, L,m), and in this section we 
comment on a few specific examples. First we take one of N, M to be zero. We start 
with the M = case, which is a natural extension of that treated in 19,. 

When m = or 1, the BAE roots assume the famous string patterns of length L—2; 
the number of strings is generically A^. We confirmed that this leads to a proper solution 
of the differential equation (|2.1j) . This means that ipM,N,mix) is, modulo a trivial change 
of variable, directly related to the ground-state expectation value of the operator Q. As 



■'"The explicit representation-theoretical construction of the Q operator is not our concern here. 
However, we mention the recent work 17 18J, where the subtleties which arise for q a root of unity are 
addressed, taking into account so-called "exact complete strings" 24 . It is worth noting that these 
strings are related to the above-mentioned possibility to multiply (I3.im by an entire function of . 

''Takemura I25| has also discussed the application of BAE to the determination of zeroes of wave 
functions for the generalised Poschl- Teller equation (or the BCi Calogero-Sutherland model), however, 
with the L2 property. The BAE itself is similar to the semi-classical form, thus different from that 
described here, the quantum form. 
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± 1.17053198009 ± 0.82128068487 i 
± 0.51014119671 ± 0.79624585390 i 
± 0.19297842512 ± 0.79284699715 i 



± 0.74741179593 ± 0.80157843231 i 
± 0.33682705394 ± 0.79395490465 i 
± 0.06295373414 ± 0.79238018877 i 



Table 1: BAE roots for L = 4, m = 0, iV = 12. 



an example, the set of Bethe ansatz roots for L = 4, m = 0, and A*" = 12 is given in 
table HI 

The case with L = 5, m = 0, A = 12 is depicted in figure ^ With increasing m, 



X XXXXXXXX X 



X XXXXXXXX X 



X XXXXXXXX X 



0.2 D.2 0.4 



Figure 1: Zeroes of Q H2.1|l at (Af, N, L, m) = (0, 12, 5, 0), which iUustrates the 3 string solution. 

the roots form longer strings and the number of roots exceeds (L — 2)N. Finally all 
but 4N zeroes are on the imaginary axis. The remaining 4N zeroes lie in 4 complex 
groups, which are empirically located near ite it i for small real part e. The BAE roots 
for {M, N, L,m) = (0,8,3,34) are plotted in figure El This behaviour will be discussed 



Figure 2: Zeroes of the wavefunction of (EU with (M, TV, i, m) (0, 8, 3, 34). 

in appendix^ in the light of an explicit solution. 

Next we consider the case M ^ and A^ = 0. When m = 0, A^ = and L odd, there 
are M (almost) strings of length L. The top roots, which are located at '^m{u) = vr, 
are displaced from the others: the distance between these roots is a little bit larger than 
others. The example {M, N, L,m) = (8,0,5,0) is shown in figure 01 This configuration 
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Figure 3: Zeroes of the wavefunction of H2.1|l with {M, N, L,m) = (8,0,5,0). We have M— 
string in the vertical direction. 

can also be interpreted as L strings of length M, rotated by 90°. The distance between 
roots in a string, however, is much less than In this interpretation, the configuration 
for the case m = 0, Nk = and L, M even is similar. For m = 0, = 0, L even and 
M odd, the pattern is a little bit different; there are L — 2 strings of length M, a string 
of length M — 1 centred at x = ni and a string of length M + 1 centred at x = 0. 

When both M and are nonzero, patterns are generally quite involved. However a 
very simple picture emerges for m = 0: the coexistence of the L — 2 strings, symmetric 
with respect to the real axis, and M strings symmetric with respect to the imaginary 
axis. The example of L = 5, M = 2, = 8, m = is shown in figure |3 



X XX^J^XX X 

ox X X^^X X X 

X xxggxx X 



Figure 4: The (M, N, L, m) = (2, 8, 5, 0) case. 



6 The field theory Umit 

It is interesting to find the field theory limit of our systems. We shall work directly 
with equation (|2.3() . and send — > oo keeping L, M and m finite. At the same time 
we focus on the region near x = by introducing a new variable z via 
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and confine our analysis to the region |z^| ^ L'^N{N + 1). Taking the large N limit, 

■4'M,N,miz) Ipmiz) with 



dz^ z dz 



z'^-'Mz) = . (6.2) 



Notice that at m = and L = 3, equation (|6.2|) coincides with the Airy equation H1.4|) . 
However, in order to identify the continuous limit of the Q function at m = 0, L = 3 
with the Airy function and hence with result of 1 from the ODE/IM correspondence, 
we should also check its asymptotic behaviour. (A possible difference in overall normal- 
isation will be ignored.) This is the point were the boundary condition ()2.5|1 becomes 
important. First notice that in terms of z the point where the condition 

ijM,N,miz)\z^zo ~ (^0 - z)^+'^ (6.3) 

was imposed (see equation (|2.5|) ) is zq = {L'^N{N + 1))^^^, so as is increased zq 
moves toward infinity. At first one might be tempted to extract the large z asymptotic 
behaviour of ipM,N,m{z) by studying the large limit of (|6.:-{|) . Since we already re- 
stricted ourselves to the region \z^\ <^ L'^N{N + 1) this would be incorrect: near the 
point z = Zq {x = \) Su linear approximation for the "potential" 

-.^.^(^^W^-^n^) 

is clearly unreasonable. However it is straightforwardly proved that, for m = and 
N > M, the condition (|6.Hj) constrains 'tpN,M,miz) to be monotonically decreasing in the 
whole segment 

< z < {L^N{N + 1))^/^ . (6.5) 

This property guarantees that the purely subdominant solution is singled out from 1)6. 2|1 . 
giving q{z) tpoiz) oc Ai(z). The argument is the following. 
The condition (|6.3|) means that 

1pM,N,m{x)\x=l-e > , tp'M,N,m{x)\x=l-e < , i^M,N,mix)\x=l~e > , (6.6) 

with a small but finite positive e. So decreasing x slightly below 1, -(/'(x) remains positive 
and in order to change the sign of ^l^'{x) the sign of •0"(x) should become negative first. 
Note now that for m = 

i^M,N,m{x) = P{x)'4}M,N,m{x) , (6.7) 

and that P{x) is positive in < x < 1 for N > M. Then the only way to have 
tp"{x) = is through ip{x) = 0. By continuity from x = 1, this contradicts the positivity 
condition ()6.6|) . 

For m > 0, due to the presence of the first derivative term in (|2.3() . this simple 
argument does not immediately apply. However, by slightly more involved reasonings 
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one can argue that at least for moderate values of m and M it is always the subdominant 
solution which is singled out in this field theory limit. For example the L = 3 and m = 1 
case related to of ^H] leads to q{x) — ^i{z) oc Ai'(x). 

Finally, we would like to mention that for m = and L general the limiting equation 
(|6.2j) coincides, up to a trivial change of variable, with the a = 1, I = 0, S = {L — 2) /2 
case of the equation 

-^Xix, E) + f (x" - Ef^ + xix, E)=0, (6.8) 

which has been identified by Lukyanov j26j with the scaling limit of the spin— XXZ 
quantum chain. 

7 Summary and discussion 

In this paper, the generalised Poschl- Teller (Heun) equation in the complex plane has 
been addressed in view of the connection relation. Remarkably, at particular values of 
parameters, a hidden link to one dimensional quantum systems of higher spins has been 
found. The Bethe ansatz method, well-developed in the theory of quantum integrable 
systems, then provides a simple characterisation of the wavefunction as an entire func- 
tion. These results begin to fill a gap in earlier studies: they show that the ODE/IM 
correspondence has a role to play in at least some finite lattice systems. The place of 
massive theories in this story, however, is yet to be clarified. 

We would like to remark that a very recent investigation of the deformed nonlinear 
a model j27j establishes a connection between perturbed Z^r parafermion theory (a 
massive theory) and the Heun equation. Though the role played by the ODE in the 
context [57j is quite far from the spirit of the ODE/IM correspondence, it would be nice 
to see whether the analysis proposed in this paper could tell us anything interesting 
about the problem ^Zl ■ 

Finally, we add some further comments on the implications of our results for the 
quantum spin chain problem. The connection rule for the ODE makes full use of the 
peculiarity of q = e*^ being a root of unity, which can be naturally extended from L = 3 
to integer values of L. Correspondingly, some peculiar features of spin model with 
S = ^ are inherited by spin models for which the quantum space possesses higher spins, 
while the spin of the auxiliary space remains at ^. Physically, vertex models, or the 
corresponding Hamiltonians, for which the quantum and the auxiliary spaces share the 
same magnitude of spin are more relevant. We call these "isotropic". Then a natural 
question arises: can we observe a similarly-simple behaviour in the largest eigenvalues 
of isotropic transfer matrices of higher spin chains, T ~ (const. )-'^'S', with a proper choice 
of coupling constant? Our preliminary numerical results answer this positively. Under 
periodic boundary conditions, the 19 vertex model (L = 4) and the 44 vertex model 
(-L = 5) show the desired simple behaviour for Ms = 3, 5, 7 when A = — cos ^. Indeed, 
for L = 5, this is confirmed by the result in section 0] and fusion relations. Although 
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these eigenvalues are characterised by the same BAE solutions as in the anisotropic 
cases, for the isotropic models they turn out to be special members of the spectrum - 
in fact, the largest in the given spin sector. We also investigated Hamiltonians with 
free-free boundaries. Through numerical diagonalisation, the spin 1 chain (L = 4) with 
quantum group invariant boundaries 28] turns out to possess the ground state energy 
Eq = — 4(A/s' — 1) for A = cos j. This is exactly the expected behaviour if T obeys the 
power law. The origin of this peculiarity, associated to the spin chain, has been argued 
for 5 = i to be the representation theory of the quantum algebra P^IHO] . There is also 
an interesting relationship between the antiferromagnetic spin-^ XXZ quantum chain 
at A = — 1/2 and a supersymmetric model of hard-core fermions . It is conceivable 
that most of the special properties emerging from our analysis will ultimately find a 
natural interpretation in the framework of similar supersymmetric systems. 

We conclude this discussion by noting that the ODE/IM correspondence has been 
extended in jn^Eniin^inSl to higher order differential equations. In these papers a 
relationship between n*'^-order ODEs and the conformal limit of SU{n) lattice models 
was established. It is interesting that these more complicated families of systems also 
possess exactly solvable points. (The corresponding differential equations are direct 
generalisations of the Airy equation (|1.4|) .) At least for these cases, the finite lattice 
extension of the models should be straightforward, and we hope to explore this point 
further in a future publication. 
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A A solution in terms of the Heun function 

The ODE ^2.1\f has an explicit solution in terms of the Heun series u = H{d, e; a, f3, 7, 6; z) 
which satisfies 

(f'u /7 5 e -.du a6(z — e) „ ,s 

with a, /3, 7, (5, e := a + /? — 7 — 5 + 1, d, e G C. (Note that e is denoted by q in M:-) 
There are four regular singularities at z = (0, 1, d, 00). It is well known that any Fuchsian 
function with four regular singular points can be transformed into the Heun function. 
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Let us make the connection between ^/J in (|B.1() and the Heun function. Set 

3;-"^/'^M,^,™(x) = —rjr^ — , {a = {m + l)/L) , (A.2) 

and adopt a variable z = = — exp(25). It is then easily established that 



i = {z-irH{d,e-a,l3,-f,6-z) , (A.3) 

with parameters 

a = a-M-N, j3 = -M - N , 7 = 1 + 0-, 5 = -2N , 

M - N 

e = -2M, d = -l, e = — — . (A.4) 

N + M ^ ' 

The case M = and m ^ 1 was numerically investigated in the main text. In this 
limit, it is immediate to see, by its degeneration to the hypergeometric function, that 

i^M,iV,m(x^) ^ (1 - X^)^. Thus 

^M,N,m{x) - , (A.5) 

which explicitly supports the asymptotic locations of the zeroes being on the unit circle, 
or on the imaginary axis in terms of u. However this does not account for the fact that 
most of them are exactly on the imaginary axis for large but finite m. 

B The determination of locations of trivial zeroes of the 
wavefunct ion 

In this appendix, we explain how to locate the 'trivial' zeroes of the wavefunction. 
Taking into account the singularities ()2.8() . there must be special solutions of the form 



'>pM,N,mix^) = , f /n , Ti-yxMi T.I') -T . ^ M,N ,m{x^) , (B.l) 



{2<y-M-N)L 

: 2 

(3;i/2 _|_ rf.-L/2^M _ 

where HM,N,miz) is nonsingular and, as a consequence of (|2.1U|) . should satisfy 

HM,N,m{ — z) OC HN,M,m{z) ■ (B.2) 

As discussed in appendix^ H]\f^M,m{z) can be written in terms of the Heun function. 
In view of this explicit form, property (|B.2() is immediate: note first that HM,N,m{x^) = 
H{—l,e,a,P,6]z). If d = —1, (|A.1|) is invariant under z — > —z, e — > — e, e <-> 5. This 
is accomplished in our case by x^ —>■ —x^, M ^ N, which can be easily verified using 
the parameterisation ()A.4|) in terms of o", M and A^. Thus the desired property ()B.2|) 
is shown to be valid. 
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Using the symmetry s — > — s of the original Poschl- Teller equation one can 

check that 

x'^^i^M^NA^/x^) (B.3) 

is also a solution of H2.3() . which turns out to be independent of HB.1|) so long as o" ^ Z. 
Therefore the particular solution which satisfies the boundary condition (|2.5j) is 

ll^M,N,m{x) = 'lpM,N,m{x^) - (-1)^ x'^^'tjjM,N,m{-^) 

1 

/ (2tT-l\l-N)L J {1\1+N)L 1 \ 

X [X 2 HM,N,mix'^) - X 2 HM,N,m{^)j ■ (B.4) 

From the expression HB.4|) one can read off the positions and orders of the trivial zeroes 
and poles. Consider first zeroes related to the symmetry (|2.11() . We denote the greatest 
common divisor of m + 1 and L by 

£ = GCD(m + l,L) , (B.5) 

and set 

IT 7T 

uj = exp(2ir/) , V = y lv' = exp(2ir/') , V = 7 ■ (B-6) 

Then at x = o;''^, {k = 0, 1, • • • ,£ — I), 'ilJM,N,m{x) has zeroes of order + 1, and at 
X = u!^, {k = 1, - ■ ■ ,L — 1) such that GCD(/c£, L) = 1, it has poles of order A^. 

There are also zeroes and poles related to the symmetries ()2.1U|) . Let ki be a positive 
integer satisfying 

{2ki + l){N + M ^)e2Z, (1<A:, <L-1). (B.7) 

By paying attention to ()B.4|) . especially the balance of the two terms in the numerator, 
we check that x = a;(2'=+i)/2 ^ pole of the order M of ipM,N,m{x) if A; / k^, while it is 
a zero of the order M + 1 when k = k^. One can easily check that if ki is a solution of 
()B.7|) then so also is /cj + j. The number of possible ki, Nf^, is thus either zero or i. 
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